\(\int \frac {1}{\sqrt [4]{a+b x^2} (c+d x^2)} \, dx\) [325]

   Optimal result
   Rubi [A] (verified)
   Mathematica [C] (warning: unable to verify)
   Maple [F]
   Fricas [F(-1)]
   Sympy [F]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 21, antiderivative size = 167 \[ \int \frac {1}{\sqrt [4]{a+b x^2} \left (c+d x^2\right )} \, dx=\frac {\sqrt [4]{a} \sqrt {-\frac {b x^2}{a}} \operatorname {EllipticPi}\left (-\frac {\sqrt {a} \sqrt {d}}{\sqrt {-b c+a d}},\arcsin \left (\frac {\sqrt [4]{a+b x^2}}{\sqrt [4]{a}}\right ),-1\right )}{\sqrt {d} \sqrt {-b c+a d} x}-\frac {\sqrt [4]{a} \sqrt {-\frac {b x^2}{a}} \operatorname {EllipticPi}\left (\frac {\sqrt {a} \sqrt {d}}{\sqrt {-b c+a d}},\arcsin \left (\frac {\sqrt [4]{a+b x^2}}{\sqrt [4]{a}}\right ),-1\right )}{\sqrt {d} \sqrt {-b c+a d} x} \]

[Out]

a^(1/4)*EllipticPi((b*x^2+a)^(1/4)/a^(1/4),-a^(1/2)*d^(1/2)/(a*d-b*c)^(1/2),I)*(-b*x^2/a)^(1/2)/x/d^(1/2)/(a*d
-b*c)^(1/2)-a^(1/4)*EllipticPi((b*x^2+a)^(1/4)/a^(1/4),a^(1/2)*d^(1/2)/(a*d-b*c)^(1/2),I)*(-b*x^2/a)^(1/2)/x/d
^(1/2)/(a*d-b*c)^(1/2)

Rubi [A] (verified)

Time = 0.08 (sec) , antiderivative size = 167, normalized size of antiderivative = 1.00, number of steps used = 4, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.143, Rules used = {408, 504, 1232} \[ \int \frac {1}{\sqrt [4]{a+b x^2} \left (c+d x^2\right )} \, dx=\frac {\sqrt [4]{a} \sqrt {-\frac {b x^2}{a}} \operatorname {EllipticPi}\left (-\frac {\sqrt {a} \sqrt {d}}{\sqrt {a d-b c}},\arcsin \left (\frac {\sqrt [4]{b x^2+a}}{\sqrt [4]{a}}\right ),-1\right )}{\sqrt {d} x \sqrt {a d-b c}}-\frac {\sqrt [4]{a} \sqrt {-\frac {b x^2}{a}} \operatorname {EllipticPi}\left (\frac {\sqrt {a} \sqrt {d}}{\sqrt {a d-b c}},\arcsin \left (\frac {\sqrt [4]{b x^2+a}}{\sqrt [4]{a}}\right ),-1\right )}{\sqrt {d} x \sqrt {a d-b c}} \]

[In]

Int[1/((a + b*x^2)^(1/4)*(c + d*x^2)),x]

[Out]

(a^(1/4)*Sqrt[-((b*x^2)/a)]*EllipticPi[-((Sqrt[a]*Sqrt[d])/Sqrt[-(b*c) + a*d]), ArcSin[(a + b*x^2)^(1/4)/a^(1/
4)], -1])/(Sqrt[d]*Sqrt[-(b*c) + a*d]*x) - (a^(1/4)*Sqrt[-((b*x^2)/a)]*EllipticPi[(Sqrt[a]*Sqrt[d])/Sqrt[-(b*c
) + a*d], ArcSin[(a + b*x^2)^(1/4)/a^(1/4)], -1])/(Sqrt[d]*Sqrt[-(b*c) + a*d]*x)

Rule 408

Int[1/(((a_) + (b_.)*(x_)^2)^(1/4)*((c_) + (d_.)*(x_)^2)), x_Symbol] :> Dist[2*(Sqrt[(-b)*(x^2/a)]/x), Subst[I
nt[x^2/(Sqrt[1 - x^4/a]*(b*c - a*d + d*x^4)), x], x, (a + b*x^2)^(1/4)], x] /; FreeQ[{a, b, c, d}, x] && NeQ[b
*c - a*d, 0]

Rule 504

Int[(x_)^2/(((a_) + (b_.)*(x_)^4)*Sqrt[(c_) + (d_.)*(x_)^4]), x_Symbol] :> With[{r = Numerator[Rt[-a/b, 2]], s
 = Denominator[Rt[-a/b, 2]]}, Dist[s/(2*b), Int[1/((r + s*x^2)*Sqrt[c + d*x^4]), x], x] - Dist[s/(2*b), Int[1/
((r - s*x^2)*Sqrt[c + d*x^4]), x], x]] /; FreeQ[{a, b, c, d}, x] && NeQ[b*c - a*d, 0]

Rule 1232

Int[1/(((d_) + (e_.)*(x_)^2)*Sqrt[(a_) + (c_.)*(x_)^4]), x_Symbol] :> With[{q = Rt[-c/a, 4]}, Simp[(1/(d*Sqrt[
a]*q))*EllipticPi[-e/(d*q^2), ArcSin[q*x], -1], x]] /; FreeQ[{a, c, d, e}, x] && NegQ[c/a] && GtQ[a, 0]

Rubi steps \begin{align*} \text {integral}& = \frac {\left (2 \sqrt {-\frac {b x^2}{a}}\right ) \text {Subst}\left (\int \frac {x^2}{\sqrt {1-\frac {x^4}{a}} \left (b c-a d+d x^4\right )} \, dx,x,\sqrt [4]{a+b x^2}\right )}{x} \\ & = -\frac {\sqrt {-\frac {b x^2}{a}} \text {Subst}\left (\int \frac {1}{\left (\sqrt {-b c+a d}-\sqrt {d} x^2\right ) \sqrt {1-\frac {x^4}{a}}} \, dx,x,\sqrt [4]{a+b x^2}\right )}{\sqrt {d} x}+\frac {\sqrt {-\frac {b x^2}{a}} \text {Subst}\left (\int \frac {1}{\left (\sqrt {-b c+a d}+\sqrt {d} x^2\right ) \sqrt {1-\frac {x^4}{a}}} \, dx,x,\sqrt [4]{a+b x^2}\right )}{\sqrt {d} x} \\ & = \frac {\sqrt [4]{a} \sqrt {-\frac {b x^2}{a}} \Pi \left (-\frac {\sqrt {a} \sqrt {d}}{\sqrt {-b c+a d}};\left .\sin ^{-1}\left (\frac {\sqrt [4]{a+b x^2}}{\sqrt [4]{a}}\right )\right |-1\right )}{\sqrt {d} \sqrt {-b c+a d} x}-\frac {\sqrt [4]{a} \sqrt {-\frac {b x^2}{a}} \Pi \left (\frac {\sqrt {a} \sqrt {d}}{\sqrt {-b c+a d}};\left .\sin ^{-1}\left (\frac {\sqrt [4]{a+b x^2}}{\sqrt [4]{a}}\right )\right |-1\right )}{\sqrt {d} \sqrt {-b c+a d} x} \\ \end{align*}

Mathematica [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 6 vs. order 4 in optimal.

Time = 7.77 (sec) , antiderivative size = 160, normalized size of antiderivative = 0.96 \[ \int \frac {1}{\sqrt [4]{a+b x^2} \left (c+d x^2\right )} \, dx=-\frac {6 a c x \operatorname {AppellF1}\left (\frac {1}{2},\frac {1}{4},1,\frac {3}{2},-\frac {b x^2}{a},-\frac {d x^2}{c}\right )}{\sqrt [4]{a+b x^2} \left (c+d x^2\right ) \left (-6 a c \operatorname {AppellF1}\left (\frac {1}{2},\frac {1}{4},1,\frac {3}{2},-\frac {b x^2}{a},-\frac {d x^2}{c}\right )+x^2 \left (4 a d \operatorname {AppellF1}\left (\frac {3}{2},\frac {1}{4},2,\frac {5}{2},-\frac {b x^2}{a},-\frac {d x^2}{c}\right )+b c \operatorname {AppellF1}\left (\frac {3}{2},\frac {5}{4},1,\frac {5}{2},-\frac {b x^2}{a},-\frac {d x^2}{c}\right )\right )\right )} \]

[In]

Integrate[1/((a + b*x^2)^(1/4)*(c + d*x^2)),x]

[Out]

(-6*a*c*x*AppellF1[1/2, 1/4, 1, 3/2, -((b*x^2)/a), -((d*x^2)/c)])/((a + b*x^2)^(1/4)*(c + d*x^2)*(-6*a*c*Appel
lF1[1/2, 1/4, 1, 3/2, -((b*x^2)/a), -((d*x^2)/c)] + x^2*(4*a*d*AppellF1[3/2, 1/4, 2, 5/2, -((b*x^2)/a), -((d*x
^2)/c)] + b*c*AppellF1[3/2, 5/4, 1, 5/2, -((b*x^2)/a), -((d*x^2)/c)])))

Maple [F]

\[\int \frac {1}{\left (b \,x^{2}+a \right )^{\frac {1}{4}} \left (d \,x^{2}+c \right )}d x\]

[In]

int(1/(b*x^2+a)^(1/4)/(d*x^2+c),x)

[Out]

int(1/(b*x^2+a)^(1/4)/(d*x^2+c),x)

Fricas [F(-1)]

Timed out. \[ \int \frac {1}{\sqrt [4]{a+b x^2} \left (c+d x^2\right )} \, dx=\text {Timed out} \]

[In]

integrate(1/(b*x^2+a)^(1/4)/(d*x^2+c),x, algorithm="fricas")

[Out]

Timed out

Sympy [F]

\[ \int \frac {1}{\sqrt [4]{a+b x^2} \left (c+d x^2\right )} \, dx=\int \frac {1}{\sqrt [4]{a + b x^{2}} \left (c + d x^{2}\right )}\, dx \]

[In]

integrate(1/(b*x**2+a)**(1/4)/(d*x**2+c),x)

[Out]

Integral(1/((a + b*x**2)**(1/4)*(c + d*x**2)), x)

Maxima [F]

\[ \int \frac {1}{\sqrt [4]{a+b x^2} \left (c+d x^2\right )} \, dx=\int { \frac {1}{{\left (b x^{2} + a\right )}^{\frac {1}{4}} {\left (d x^{2} + c\right )}} \,d x } \]

[In]

integrate(1/(b*x^2+a)^(1/4)/(d*x^2+c),x, algorithm="maxima")

[Out]

integrate(1/((b*x^2 + a)^(1/4)*(d*x^2 + c)), x)

Giac [F]

\[ \int \frac {1}{\sqrt [4]{a+b x^2} \left (c+d x^2\right )} \, dx=\int { \frac {1}{{\left (b x^{2} + a\right )}^{\frac {1}{4}} {\left (d x^{2} + c\right )}} \,d x } \]

[In]

integrate(1/(b*x^2+a)^(1/4)/(d*x^2+c),x, algorithm="giac")

[Out]

integrate(1/((b*x^2 + a)^(1/4)*(d*x^2 + c)), x)

Mupad [F(-1)]

Timed out. \[ \int \frac {1}{\sqrt [4]{a+b x^2} \left (c+d x^2\right )} \, dx=\int \frac {1}{{\left (b\,x^2+a\right )}^{1/4}\,\left (d\,x^2+c\right )} \,d x \]

[In]

int(1/((a + b*x^2)^(1/4)*(c + d*x^2)),x)

[Out]

int(1/((a + b*x^2)^(1/4)*(c + d*x^2)), x)